 
%WinBUGS wavelets
clear all
close all


  disp('USO Example')
  lw = 2; 
  set(0, 'DefaultAxesFontSize', 16);
  fs = 15;
  msize = 5;
  

%-------------------------------
    x1=0:0.02:0.7;
    r1=(1-betacdf(x1,2,4))./(1-betacdf(x1,3,6));
  figure(1)
  plot(x1,r1, 'linewidth', lw  )
  axis([0 0.7 0.9 1.3])
  print -deps 'C:\NPBook\ProbBasics\USOgraph.eps'
  
    x2=0.01:0.01:.99;
    r2=(betapdf(x2,2,4))./(betapdf(x2,3,6));
  figure(2)
  plot(x2,r2, 'linewidth', lw )
   axis([0 1 0 70])
  print -deps 'C:\NPBook\ProbBasics\LROgraph.eps'